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COMMUNICATION CHANNEL OPTIMIZATION SYSTEMS AND METHODS IN MULTI- 
USER COMMUNICATION SYSTEMS 

Cross -Reference to Related Applications 

This application claims priority from United States 
5 provisional patent applications Serial Nos. 60/517,389, filed 
on November 6, 2003, and 60/517,893, filed on November 7, 2003. 
The entire contents of each of these provisional applications 
are hereby incorporated herein by reference. 

Field of the Invention 

10 The invention relates to communication systems in 

general, and particularly to communication channel optimization 
in multi-user communication systems. 

Background 

In any communication system, the quality and capacity 
15 of a communication channel are affected by such factors as 
interference, allocation of communication resources, the 
communication schemes or algorithms used on the communication 
channel, and the particular communication equipment implemented 
at transmitting and receiving ends of the channel. 

20 The effects of certain factors may be mitigated 

through efficient resource allocation and selection of 
communication schemes and equipment. According to some 
conventional communication techniques, processing operations 
intended to compensate for other communication channel effects 

25 are primarily receiver based. For example, interference 

cancellation is performed by a receiver in known communication 
systems. In addition, the implementation of different types of 
communication equipment in conjunction with the same type of 
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channel, such as different communication terminals in a 
wireless communication system for instance, may affect received 
signal processing operations at all receivers. 

Communication terminals at the ends of a 
5 communication channel are seldom identical. In wireless 
communication systems, for example, user communication 
terminals at one end of a communication channel normally have 
much more limited resources than base stations. In known MIMO 
(Multiple Input Multiple Output) systems, each receiver has at 

10 least as many antennas as a transmitter. This constraint is 
difficult to satisfy where communication equipment on opposite 
ends of a communication channel are significantly different,, as 
in the case of wireless communication terminals and base 
stations in wireless communication systems, for example. In 

15 addition, resource limitations at one receiver in such a multi- 
user system can also affect other receivers in the system. 

Summary of the Invention 

According to one aspect of the invention, a method of 
processing signals to be transmitted to receivers on 

20 communication channels is provided. The method includes 

determining pre-coding signal weights based on channel state 
information associated with the communication channels to 
provide proportional power allocation to the signals, and 
applying the signal weights to the signals. In a preferred 

25 embodiment, the channel state information is received from the 
receivers . 

The invention also provides, in another aspect, a 
method which includes receiving over a sub-group of 
communication channels a subset of signals to which pre-coding 
30 signal weights based on channel state information associated 
with the communication channels to provide proportional power 
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allocation have been applied. The received subset of signals 
using inverses of the pre-coding signal weights based on 
channel state information associated with the sub-group of 
channels to decode the received subset of signals. 

According to an embodiment of the invention, the 
signals to be transmitted include respective groups of signals 
to be transmitted to the receivers, and the pre-coding weights 
are determined to separate the respective groups of signals. 
Decoding then separates individual signals in the received 
subset of signals. 

In another aspect, the invention provides a system 
for processing signals to be transmitted to receivers on 
communication channels. The system preferably includes an 
input for receiving the signals and a processor. The processor 
is configured to determine pre-coding signal weights based on 
channel state information associated with the communication 
channels to provide proportional power allocation to the 
signals, and to apply the signal weights to the signals. In a 
preferred embodiment, the system further includes multiple 
antennas which provide respective sub-MIMO channels to the 
receivers . 

The invention also provides a system that includes an 
input an input for receiving over a sub-group of communication 
channels a subset of signals to which pre-coding signal weights 
25 based on channel state information associated with the 

communication channels to provide proportional power allocation 
have been applied, and a processor. The processor is 
configured to decode the received subset of signals using 
inverses of the pre-coding signal weights based on channel 
30 state information associated with the sub-group of the 
channels . 
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A further method of processing signals to be 
concurrently transmitted to receivers over communication 
channels, in accordance with still another aspect of the 
invention, includes determining channel state information for 
5 the communication channels, determining a spatial coding matrix 
which includes respective set of spatial coding weights for 
each of the receivers based on the channel state information, 
and applying the spatial coding weights in the spatial coding 
matrix to the signals . 

10 A method in accordance with a still further aspect of 

the invention includes determining channel state information 
for a communication channel between a receiver and a 
transmitter, transmitting the channel state information to the 
transmitter, and receiving from the transmitter one of a 

15 plurality of demodulation matrices for demodulating 

subsequently received communication signals to which spatial 
coding weights comprising respective sets of spatial coding 
weights for multiple receivers have been applied. 

A network element for processing signals to be 
20 concurrently transmitted to multiple communication terminals in 
a communication network is also provided. The network element 
preferably includes an input configured to receive the signals, 
and a processor. The processor is configured to determine 
channel state information for each communication channel 
25 between the network element and the communication terminals, to 
determine a spatial coding matrix comprising a respective set 
of spatial coding weights for each of the communication 
terminals based on the channel state information, and to apply 
the spatial coding weights in the spatial coding matrix to the 
30 signals. 

In a related aspect, a communication terminal for 
operation in a communication network is provided. A processor 
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in the terminal is configured to determine channel state 
information for communication channels between the 
communication terminal and a network element in the 
communication network. The terminal also includes at least one 
5 antenna for transmitting the channel state information from the 
communication terminal to the network element, receiving a 
demodulation matrix from the network element, and receiving 
signals concurrently transmitted to multiple communication 
terminals from the network element. The processor is further 
10 configured to demodulate the received signals using the 
demodulation matrix . 

Other aspects and features of the present invention 
will become apparent to those ordinarily skilled in the art 
upon review of the following description of the specific 
15 embodiments of the invention. 

Brief Description of the Drawings 

The invention will now be described in greater detail 
with reference to the accompanying diagrams, in which: 

Fig. 1 is a block diagram of a 6x6 MIMO system; 

20 Fig. 2 is a block diagram of a decomposed 6x6 MIMO 

system; 

Fig. 3 is a block diagram of a MIMO system, 
illustrating inter-user interference; 

Fig. 4 is a block diagram of a MIMO system using 
25 polarized communication channels, illustrating inter-user 
interference; 

Fig. 5 shows a block diagram of a MIMO system 
according to an embodiment of the invention; 



Fig. 6 is a block diagram of a MIMO BLAST system; 
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Fig. 7 is a block diagram of a multi-user MIMO system 
in accordance with another embodiment of the invention; 

Fig. 8 is a block diagram of a conventional null 
beamforming MIMO system; 

5 Fig. 9 is a block diagram of a known multi-user BLAST 

MIMO system; 

Fig. 10 is a block diagram of a known round-robin TDM 
(Time Division Multiplexing) BLAST MIMO system; 

Fig. 11 is a plot of BLER (Block Error Rate) versus 
10 SNR (Signal-to-Noise Ratio) for an embodiment of the invention 
and several known communication schemes; 

Fig. 12 is a plot of BLER versus Eb/No (Energy per 
Bit to Spectral Noise Density ratio) for an embodiment of the 
invention and several known communication schemes; and 

15 Fig. 13 is a plot of BLER versus SNR for several 

embodiments of the invention. 

Detailed Description of the Preferred Embodiments 

According to embodiments of the present invention, 
systems and methods are provided which enhance the performance 
20 of communication channels in a communication system, to thereby 
improve, for example, the transmission performance of multi- 
user MIMO (Multiple Input Multiple Output) communication 
systems . 

In MIMO systems, a multi-data stream transmitter at 
25 a base transceiver station (BTS) that provides communication 
services for a coverage area or cell in a wireless 
communication system transmits communication signals to user 
terminals via multiple antennas. User terminals are also 
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commonly referred to as user equipment (UE) , communication 
devices, and mobile stations, for instance. At a UE side, 
multiple receive antennas are employed for each user. 

Fig. 1 is a block diagram of a 6x6 MIMO system, 
5 which includes a plurality of antennas 16 at a BTS 10, and a 
plurality of antennas 18 and a MIMO decoder 20 at a UE 12. In 
this system, 6 communication signals 14 intended for the UE 20, 
labelled as Si (1) through s 6 (1) are transmitted via the antennas 
16 from the BTS 10 to the UE 12. At the UE 12, each of the 
10 antennas 18 receives the signals transmitted from the antennas 
16, and the received signals are decoded in a MIMO decoder 20. 

It should be appreciated that the system of Fig. 1 
is intended for illustrative purposes only. As will be 
apparent to those skilled in the art to which the present 

15 invention pertains, BTSs include further components in addition 
to the antennas 16, such as components to generate the signals 
si (1> through s 6 (1) for instance. Similarly, the UE 12 includes 
components to further process received signals decoded by the 
MIMO decoder 20. Also, the BTS 10 and the UE 12 normally 

20 support both transmit and receive operations. 

Known MIMO systems do not support simultaneous multi- 
user MIMO transmissions where each UE does not have at least 
the same number of antennas as a BTS. Instead, MIMO is 
typically used as a single "fat-pipe", and multiple users are 

25 served through the use of time division techniques. In 

addition, it is practically difficult to realize very large 
dimension MIMO systems, 8x8 systems for example, due the 
physical size limitations of UEs. As the physical size of a UE 
is usually limited, the distance by which multiple antennas for 

30 larger dimension MIMO systems can be separated is also limited, 
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such that the antennas at the UE become highly correlated. 
This drastically reduces the MIMO channel capacity. 

However, large dimension MIMO systems may be 
decomposed into combined sub-MIMO systems, because in general, 
5 channel fading between different UEs is un-correlated. MIMO 
channel capacity can then be more efficiently exploited. Fig. 
2 is a block diagram of a decomposed 6x6 MIMO system. The 
decomposed 6x6 MIMO system includes antennas 32 at a BTS 22, 
and three UEs 24, 26, 28, each of which includes a pair of 

10 antennas 34/36, 40/42, 46/48 and a MIMO decoder 38, 44, 50. At 
the BTS 22, signals Si (1) /s 2 (1) , Si <2) /s 2 <2) , Si (3) /s 2 (3) are 
transmitted via a respective antenna in each of a plurality of 
sub-groups, pairs in Fig. 2, of the antennas 32 to 
corresponding UEs 24, 26, 28. The superscripts (1), (2), (3) 

15 indicate for which UE 24, 26, 28, a signal is intended. Each 
pair of the antennas 32 forms a 2x2 sub-MIMO channel with the 
pair of antennas of a respective one of the UEs 24, 26, 28. 

To apply the MIMO technique to a multi-user system, 
inter-user MIMO interference is a major issue. In Fig. 2, for 

20 example, although each pair of signals^ Si (1) /s 2 (1) , Si (2) /s 2 (2> , 
si (3) /s 2 <3) may be intended for a particular UE 24, 26, 28, the 
antennas 34/36, 40/42, 46/48 in each of the UEs 24, 26, 28 
receive communication signals from all of the antennas 32. In 
order to cancel such interference using known techniques, the 

25 number of receive antennas must be equal to or greater than the 
number of transmit antennas. Due primarily to physical space 
and form factor constraints on UEs, the number of antennas that 
can be provided at a UE is limited. Thus, for downlink (BTS to 
UE) transmissions, the number of receive antennas is typically 

30 smaller than the number of transmit antennas. 
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Figs. 3 and 4 are block diagrams of MIMO systems, 
and illustrate inter-user interference. In the system of Fig. 
3, the signals 58, labelled Si a> /s 2 (1) , Si (2) /s 2 (2) , are transmitted 
from the BTS 52 via respective ones of each pair of the 
5 antennas 60 to the UEs 54, 56. Signals received by the 

antennas 62/64, 68/70 in the UEs 54, 56 are processed by the 
MIMO decoders 66, 72. Interference in the MIMO system of Fig. 
3 is indicated at 73. As shown, any communication signals that 
are received at one of the UEs 54, 56 but intended for the 
10 other of the UEs 54, 56 represent interference at that UE . For 
example, versions of si <2) and s 2 <2) received at the UE 54 
represent interference . 

Similarly, in Fig. 4, the signals 80, labelled 
si (1) /s 2 (1) , si (2) /s 2 <2) , are transmitted from the BTS 74 via 

15 respective ones of each pair of the antennas 82 to the UEs 76, 
78, received by the antennas 84/86, 90/92, and processed by the 
MIMO decoders 88, 94. Polarized channels, provided by the 
vertical- and horizontal- polarization antennas V and H in the 
example system of Fig. 4, reduce interference 95 relative to 

20 the general MIMO system of Fig. 3, but do not cancel the inter- 
user interference to such a degree that interference 
cancellation is unnecessary. 

The task of interference cancellation is typically 
performed at a user terminal. In accordance with an embodiment 

25 of the invention, downlink communication channel interference 
cancellation is effectively split between the BTS (transmit) 
side and the UE (receive) side. For example, a BTS may perform 
inter-user separation based pre-coding of data to be 
transmitted, while at the UE side, a UE performs MIMO layer 

30 separation and decoding. 
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One type of multi-user MIMO system in which the 
invention may be implemented delivers communication signals 
according to the layered space-time known as MIMO-BLAST 
concurrently to multiple users. Such a system is preferably 
5 realized with feedback of communication channel state 

information from each UE to a BTS . Channel state feedback 
techniques are very well suited for application in conjunction 
with fixed or nomadic wireless communication channels due to 
the slow variation of such channels, which allows accurate 
10 channel state information feedback from the UEs to the BTS. 

In one embodiment, a closed-loop pre-coded transmit 
antenna array for sub-MIMO transmission, preferably MIMO-BLAST 
transmission, in a multi-user environment is provided. Channel 
state information is measured by or fed back to a BTS, and at 

15 the BTS side, jointly optimized weights are computed and 
applied to antenna input signals to cancel inter-user MIMO 
interference. Therefore, in one sense, a system in accordance 
with this embodiment of the invention may be considered as an 
adaptive weighted transmit antenna array operating in the MIMO- 

20 BLAST mode. This concept is a departure from the conventional 
beamforming phased antenna array. 

A MIMO system can be expressed as 

y = Hs+f], (1) 

where 

25 y = [y, y 2 ... y N ] T is a vector of communication signals 

received at a receiver; 



s = i s \ s i s m Y ^ s a vector of communication signals 

transmitted by a transmitter; 
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V = ^2 InY ^ s a vector of noise components affecting 

the transmitted communication signals; 



h 2V h 22 



L 2M 



NM 



is a channel matrix of communication 



channel attenuation factors; 

N is a number of antennas at the receiver; and 

M is a number of antennas at the transmitter. 

To decode the transmitted signal s , the receiver 
performs the inverse process 



s^Gy-Gf) , 



(2) 



10 where 



G = H + =(H'H)- l H' 



(3) 



is the Moore-Penrose pseudo-inverse of H , and H' is a 
conjugate matrix of H , illustratively a Hermitian conjugation 
or conjugate transpose. The post-detection SNR (signal-to- 
15 noise ratio) for a decoded element s f . of s is given by 



" , ' = - 2 tl^ll 2 ' 



(4) 



where g, is the / th row of G , and a 2 represents the variance 
of the elements in if (assuming that all the elements in fj have 
the same level of variance) . From equation (4) , it can be seen 
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that the SNR of post-detection signals is determined by G , 
which is given by equation (3) . Note that ||g,.|| 2 is not only 
determined by channel attenuation factors |/*, y | , but by the 
condition of the channel matrix H as well. When H is ill- 
5 conditioned, can be very large; hence y ( will be very 

small . 

From equation (2), it can be seen that the post- 
detection signal power is a fixed value |s f .| 2 , such that y i is in 
fact determined by the second term Gfj only. System performance 
10 can be improved by reducing post-detection noise, which is 
represented by ||g.|| 2 . One possible issue to consider in this 

regard is whether, when a transmitter has information about H , 
pre-equalization can improve system performance. 

By defining a new signal vector, illustratively for 
15 the square channel matrix case in which M = N , 

x = Js=H r {HH r Ys f (5) 
we have 

r=Hx = s. (6) 

The matrices G and J are equal in this example, 
20 when H is a square matrix (M = N). 



The SNR for this pre-equalized signal can be 
expressed as 
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*/=riT — ' (7) 



where j t is the ith row of J . 

The pre-equalization approach is similar to power 
control, i.e., the weak user gets more power so that all the 
5 users are equal. However, this is not an efficient approach to 
utilize system power. 

Another important observation is that the pre- 
equalization matrix J is completely determined by the channel 
matrix H . That is, except for making up an identity matrix, 
10 there are no other kinds of optimizations in J . 

An embodiment of the invention provides for user 
separation at the transmitter. One preferred user separation 
technique allows the use of ML (Maximum Likelihood) detection 
schemes at a receiver, such that the diversity order for each 
15 receiver is increased. In addition, since layers need not be 
separated within a BTS, system performance may be improved in 
at least two further aspects, namely, to enhance an equivalent 
channel matrix and proportional transmitting power allocation. 

Fig. 5 shows a block diagram of a MIMO system 
20 according to an embodiment of the invention. The system 

includes a BTS 100 having a pre-coder 106 and a plurality of 
antennas 108, 110, 112, 114, and UEs 102, 104, each having a 
plurality of antennas 116/118, 122/124 and a MIMO decoder 120, 
126. In Fig. 5, M = 4 , AT . = 2 is the number of antennas at an 
25 ith UE, U = 2 is the number of UEs, and M = U*N i . In a 

particularly preferred embodiment, the system of Fig. 5 is a 
MIMO BLAST system. 
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Of course, the system of Fig. 5 is one illustrative 
example of a system in which the invention may be implemented. 
The invention is in no way limited thereto. Extension of the 
principles of the present invention to systems having other 
5 dimensions will be apparent to those skilled in the art. 

At the BTS side, the BTS 100 preferably uses the U 
degrees of freedom (one per UE 102, 104) of the transmit 
antennas 108, 110, 112, 114 to perform weighted pre-coding of 
the signals Si (1) , s 2 (1) , Si (2) , S2 (2> in the pre-coder 106, while 
10 reserving the N. degrees of freedom (one per receive antenna of 
each UE 102, 104) of the transmit antennas 108, 110, 112, 114 
to maximize the N i xN i sub-MIMO channel capacity and 
proportional antenna power allocation for each user. 

At the UE side, each UE 102, 104 decodes 
15 communication signals received at its antennas 116/118, 
122/124, using ML or MMSE (Minimum Mean Squared Error) 
decoding, for example. The N i antennas also provide diversity 
gain . 

The UE-1 102 preferably determines and feeds back 
20 channel state information H x (4x2) to the BTS 100. The UE-2 104 
similarly preferably determines and feeds back channel state 
information H 2 (4x2) to the BTS 100. For some types of 
communication channel, channel state information may instead be 
determined locally by the BTS 100. Thus, channel state 
25 information determination is shown conceptually in Fig. 5 
outside the UEs 102, 104. Those skilled in the art will 
appreciate that channel state determination may be performed in 
the UEs 102, 104 or the BTS 100 by a digital signal processor 
(DSP) or a general-purpose processor adapted to execute signal 
30 processing software, for example. Various techniques for 

determining channel state information will be apparent to those 
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skilled in the art, including estimation based on pilot 
channels in CDMA (Code Division Multiple Access) systems or 
preambles and scattered pilot tones in OFDM (Orthogonal 
Frequency Division Multiplexing) systems, for example. 

5 Based on the channel state information Hi and H 2 , the 

BTS 100 computes an antenna weight matrix or pre-coding matrix 
P t which preferably cancels inter-user MIMO interference 
between UE-1 102 and UE-2 104 and maximizes MIMO system channel 
capacity for UE-1 102 and UE-2 104. 

10 Before proceeding with a detailed analysis of the 

system of Fig. 5, a combined (with respect to users) 2x2 MIMO 
system, in which M = 2, N t = 1 , and U = 2 , is first considered. 
In this system, there are no multiple layers associated with 
each user. The channel matrix H can be expressed as 



15 



H = 



'Ml "12 
_^21 ^22. 



(8) 



Now we define a pre-coding matrix P , so that 



C = HP = 



K KiXpix P22 



h uPu + KlPlK h \\P\2+ h nP22 

h 2i p u +h 22 p 2i h 2] p l2 +h 22 p 22 



(9) 



One goal is to identify if a solution to the following eguation 
exists : 



20 



{max\h liPu +h n p 2l \ 

< P\\*Pz\ 

[h 2x p u +h 22 p 2l =0 
max\h 2j p y2 +h 22 p 22 \ 
h nP i 2 +h nP22 =° 



(10) 
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In the pre-equalization approach, the set of elements 
in a pre-equalization matrix £/, i » 721 » V 12 » 7*22 } are used to satisfy 
the following condition: 

[ h \J\\ ^KiJix = 1 
\ h 2\J\ \ + h 2ih\ =° 



5 Given the same ||j3 f || 2 = |/ ( | 2 / with 1 = 1,2, equation (11) forces 
h n j xx +h xl j 2X and h lx j xl + h 21 j 22 to 1, while equation (10) tries to 
maximize the power of the analogous components h u p u + /*, 2 P 2 i and 
^2\P\2 + ^22^22 • This illustrates one primary difference between 
the user separation techniques according to embodiments of the 
10 invention and pre-equalization approaches. 

Equation (10) can be manipulated to the form of 



max 

Pu 



M2I 



'22 y 



P21 =~ 



'22 



max 

P22 



^22 



h 22 - 



h 2l h l2 



'11 y 



(12) 



h l2 p 22 



Although there is no optimized solution to equation 
(12) , selection of p u and p 22 such that 



15 



Pu 



oc 



'22 



and |p 22 |oc 



(13) 
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improves the system capacity. 

If each UE has only one receive antenna and receives 
only one layer of the MIMO signal, then proportional power 
allocation may be achieved. However, if a UE has multiple 
5 antennas and receives multiple layers of signals, as in Fig. 5, 
then the situation will be different. 

As described briefly above, the pre-coder 106 at the 
BTS 100 determines and applies pre-coding weights to the 
signals Si (1) , s 2 (1) , si (2) , s 2 <2) to perform the function of inter- 

10 user interference cancellation through user separation. Each 

receiver, UEs 102, 104 in Fig. 5, then performs the function of 
inter-antenna interference cancellation. Thus, algorithms such 
as ML and iterative ZF (Zero Forcing) /MMSE can be used to 
further improve detection results. The concept of splitting 

15 interference cancellation functions between the transmitter and 
receivers is described in further detail below. 

For the combined 2x2 MIMO multi-user system of Fig. 
5, the pre-coding matrix P is selected so that the combined 
channel matrix C has the format of 





C,2 


0 


0 


c 2 . 




0 


0 


0 


0 


^33 


C 34 


0 


0 


^43 


C 44 



(14) 



Note that only c, and c 2 , where c t represents the i-th row of C, 
affect the UE-1 102, and that only c 3 and c 4 affect the UE-2 
104. The combined channel matrix C has a form of U N i ^N i 

sub-matrices, diagonal elements of which are respective 
25 diagonal elements of C . Elements of C outside the plurality 
of N.xNj sub-matrices are zero. In another sense, C may be 
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considered as having groups of rows associated with UEs and 
groups of columns respectively associated with antenna pairs. 
For example, c, and c 2 are associated with the UE-1 102 as 
described above, and the columns of C are respectively 
5 associated with the antennas 108, 110, 112, 114. As the 

antenna pairs 108/110 and 112/114 provide sub-MIMO channels to 
the UEs 102 and 104, each pair may be considered to be 
associated with a respective one of the UEs 102 and 104. Thus, 
each element of C positioned in a row associated with a 
10 particular UE and a column corresponding to an antenna that is 
associated with a different UE is forced to zero by selection 
of pre-coding weights in the pre-coding matrix P . 

The first two columns of the pre-coding matrix P are 
determined such that they satisfy 



15 



C M =K\P\\ + h \lPl\ + / *13/>3. +KaP4\ 
C 2\ = h 2\P\\ + h 2lPl\ + h 2)Pl\ + h 2AP*\ 
KP\\ +KP2X + / *33/ ? 3l +KPAX =0 
KPu +^42^21+^43^31 + KPa\ =0 



, and 



(15) 



<?i2 =h uPn +h l2 p 22 +h l3 p n +h X4 p 42 

^ C 22 = K\P\2 + ^22^22 + ^23^32 + KaPa2 
KP\2 + ^32^22 + ^33^32 + KPa2 = 0 
KPm + k ^P22 + KPz2 + KPA2 = 0 



The second two columns of P, related to the UE-2 104, are 
preferably determined in an analogous manner. 

From equation (14), it is not difficult to see that 
20 the equivalent system for UE-1 102 is 









y 2 . 




_ C 21 C 22JL' S 2. 



(17) 
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which represents a 2x2 MIMO system. When an ML decoder is used 
at a UE, the diversity order is two. 

The decoders 120, 126 decode received signals using 
an inverse matrix such as the Moore-Penrose pseudo-inverse 
5 matrix of a corresponding sub-matrix of P . In Fig. 5, Q } and 
Q 2 indicate sub-matrices of such an inverse matrix Q that 
relate to the UEs 102, 104, respectively. Preferably, an 
inverse matrix D of C = HP , or strictly sub-matrices D x and D 2 
thereof, are used by the decoders 120, 126 in Fig. 5. These 
10 matrices can be derived by each receiver based on channel 
estimation using pilot channels, for . example. In one 
embodiment, the BTS 100 sends pilot tones to UEs 102, 104, and 
each UE the estimates the elements of its corresponding inverse 
sub-matrix . 

15 For brevity, the following analysis relates only to 

the UE-1 102. Those skilled in the art will appreciate that 
the analysis for the UE-2 104 would be similar. 

In equation (15), let p 3l and p 41 force c 13 = c, 4 = 0 , 
which gives 



20 





Pu 


1 






A 


where 


A = 





-h 



43 



*34'*43 



"32 
41 Kl 



Pu 

IP2\} 



(18) 



If A is defined as 



A = 



h 2X h 22 



"44 
h 



43 



'34 



*33 



^3. h 32~\ 
K\ A 42_T 



then (19) 
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= A 


Pn 


- C 2\_ 




.P21. 



(20) 



As the matrix A is determined by the channel matrix H , p u and 
p l2 are free to be chosen, and as such can be used for channel 
matrix optimization . 



P32 

P*2 J 



Similarly, from equation (16), 

and 



1 






A 






= A 




_ C 22_ 




.P22. 



A 34~p3. ^32 
'33 1*41 ^42. 



P\2 
.P22J 



(21) 



(22) 



By combining equations (20) and (22) , we can 



establish a relation between 





c 12 " 


and 


Pn 


Pn' 






_ C 21 


C 22. 




-P21 


P22. 



as follows: 



10 





= A 


Pn 


P\2 


_ C 2\ C 22. 




P2\ 


P22. 



(23) 



and the parameters in 



PU P\2 
Pl\ Pl2_ 



channel matrix 



L C 21 C 22j 



case 



[:;] 



provide for optimization of the 
Note that in the pre-equalization 
is selected in such a way that 
On the contrary, in this procedure, two goals are 



Pl\ P\2 
P2\ P22 



C U C \2 
L C 2I C 22_ 



is set to 



15 achieved during pre-coding, namely, separating the layers with 
respect to receivers and allocating transmitting power to 
facilitate equal layer performance. Since individual layers 
now no longer need to be separated at a transmitter, and power 
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allocation is not pre-equalization, the matrix 



P22J 



may be 



chosen according to different criteria, such as improving 
channel matrix condition and providing proportional power 
allocation, for example. 

To improve channel matrix condition, 

should be maximized. One possible way to achieve this is to 
force the elements in A to add constructively to form the 
diagonal elements c n and c 22 . In particular, the pre-coding 
weights for the UE-1 102 may be set to 





f 




\ 


det 


C \l C \2 






V 


_?2\ C 22_ 


J 



10 



P21 = va n 



[P22 = ™ 22 



(24) 



where v is a power normalization factor and the elements a 0 are 
elements of A. Substituting (24) into (23) yields 



+ \a 



= v 



a u a 2l +a n a 22 



2l"ll +a 22 a i2 



4 21 



*22 



(25) 



It can thus be seen that c u and c 22 are enhanced 
15 constructively, while c X2 and c 21 are constructed randomly. 
Therefore, the condition of C becomes more robust. From a 
beam-forming point of view, layer-1 Si (1) and layer-2 s 2 (1) are 
beamed onto antenna-1 108 and antenna-2 110, respectively, 
following the MRC (Maximum Ratio Combining) criterion. The 
20 elements c, 2 and c 2 , represent both inter-layer interference and 
receiver diversity. Recall that X x + X 2 = c u + c 22 , where X. (i = 1,2) 
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are eigenvalues of ° u ° n . Since c u and c 22 are enhanced 

L C 21 c n\ 

constructively, so are A,+/l 2 . In addition, according to 



Cauchy-Schwarz Inequality, we always have 



a l2 * a 22 





r 




\ 




~c u c l2 ~ 


det 








V 


_ C 2\ C 22_ 


) 



> 0 when 



5 For the UE-2 104, the elements of P may be selected 

in an analogous manner. The elements p u , /? 23 , p l4 , and p 24 are 
preferably selected to force c 3l = c n = c 4l - c 42 =0 , and /? 33 , p 43 f p 34 , 
and are preferably selected as P33 = va 33 , p 43 = va\ 4 , p M =va* 43 , 

and p 44 =ya* 44f where v is as defined above. For UE-2 104, 



10 however, A = 



A = h u h 22 -h l2 h 2l . 



"A33 


K~ 


1 




/j 32 " 




-h X2 










A 


A. 


^42. 


-h 2X 


K . 




h 24_ 



and 



A new scheme to further enhance multi-user MIMO 
system performance has been described. The above embodiments 
are based on splitting the interference cancellation task 

15 between a transmitter and receivers. Specifically, a 

transmitter performs inter-user separation pre-coding to define 
sub-MIMO channels, while receivers perform individual layer 
separation. The transmitter and each receiver benefit from 
this task partitioning. From the transmitter point of view, 

20 since it is no longer required to provide individual layer 

separation, it has the freedom to perform beamforming, which 
results in more robust equivalent channel matrix, and 
proportional power allocation. At the receiver, since multiple 
antennas are receiving signals from multiple layers, when the 

25 ML decoding algorithm is used, additional diversity gain can be 
achieved. 
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The above description relates primarily to MIMO 
systems in which sub-MIMO channels between the BTS 100 and each 
UE 102, 104 have the same dimension. However, it should be 
appreciated that these embodiments of the invention may also be 
5 extended to systems in which UEs do not have the same number of 
antennas, such that sub-MIMO channels of different dimensions 
are supported in the same system. 

Further embodiments of the present invention will be 
best appreciated in conjunction with the following detailed 
10 analysis of MIMO, particularly MIMO BLAST. Fig. 6 is a block 
diagram of a MIMO BLAST system. 

The system of Fig. 6 includes a BTS 130 and a UE 132. 
The BTS includes a modulator 134, an S/P (serial-to-parallel) 
converter 136, and a plurality of M antennas 140. At the UK 
15 132, communication signals received at a plurality of TV 
antennas 142 are decoded in a decoder 150 to recover the 
transmitted signals si (1) ,. ., s M (1) at 152. Communication 
channel noise is represented by the adders 144, 146, 148. 

Fig. 6 shows a general case of an MxN MIMO system. 

20 Although three representative signals and antennas are shown at 
the BTS 130 and the UE 132, the invention is in no way 
restricted any particular dimension of MIMO system. M and N 
may be equal to, larger than, or smaller than 3. Also, as 
described above in conjunction with Fig. 1, BTSs and UEs 

25 typically include further components that have not been shown 
in Fig. 6 to avoid congestion in the drawing. 

According to the basic point-to-point Layered STC 
architecture known as BLAST, a sequence of the modulated 
symbols, each symbol having a duration of T , from the modulator 
30 134, is serial-to-parallel converted in the S/P converter 136 
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into parallel transmission blocks in the signals 138. Each 
transmission block consists of K ch symbols, where K ch is the 
number of spatial channels. In the downlink transmission case, 
K ch is equal to the number of antennas, M , at the BTS 130. 
5 All of the symbols of a transmission block are simultaneously 
radiated into space, and each symbol is radiated by a 
respective one of the antennas 140. As the equivalent time 
duration of the transmission block relative to the original 
modulated signal output by the modulator 134 is K ch T , the 
10 radiated signal requires spectrum width equal to a fraction 
(\l K ch ) of that of the original signal. This achieves a very 
high spectral efficiency. 

The reception of a signal in such systems requires 
multiple antennas, with N>M. The model of a received signal 
15 is described by following vector-matrix expression 



y = ^PJMHs+fj , (26) 

where 

P s is the total transmitted power; and 

y, M , H , s, and fj are as defined above. 

20 The elements of H are preferably independent 

Gaussian complex random variables, with zero mean and isf)/*^! 2 } = 1 
variances. fj also preferably has zero mean, and R-2<j^I N 
covariance matrix, where I N is an NxN -dimensioned unit matrix 
and a 2 is the variance of one quadrature component of fj . The 
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mean power of each radiated symbol is equal to unity, i.e., 
E<\s m \ 2 } = l. 

The solution for linear estimation of a vector of 
modulated symbols can be carried out, for example, by the ZF 
5 criterion expressed as 

S =H+y/(P s /M) . (27) 

Estimation by the MMSE criterion is as follows: 

S =(H'H/N + (2a n 2 1 P s )I M )' x {H'4M)y 

= (H'jM)(HH'/N + (2a n 2 /P s )I M y x y , (28) 



10 where I M is an M x M -dimensioned unit matrix 



The MMSE algorithm provides a significant gain as 
contrasted with the ZF algorithm in a channel with Rayleigh 
fading . 



For single user point-to-point MIMO with an open-loop 
15 transmission, the MIMO channel capacity is proportional to 

min{N,M} . It well known that the throughput capacity of MIMO 
grows linearly with an increase of min{Af,M} . Conventional 
MIMO-BLAST receiver processing schemes require co-processing of 
signals received by all antennas in the MIMO system. However, 
20 such a method cannot be applied for multi-user combined MIMO 
systems, as the signals received by all other UEs may not 
always be accessible to each UE. Therefore, usage of MIMO 
BLAST in multi-user systems to provide for multiple access may 
be inefficient, despite a marginal increase of throughput 
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capacity at the expense of an increase in the number 
transmitting and receiving antennas. 

In the multi-user point-to-multi-point case, since 
the inter-user communication is typically not done, the 
5 capacity of a multi-user system with open-loop transmit 

diversity is determined by min{M, N { , N 2 >...N U } , where U , as above, 
is the number of UEs in a multi-user system. Thus, common 
throughput capacity will be limited by UEs with the least 
number of receiving antennas. 

10 Fig. 7 is a block diagram of a multi-user MIMO system 

in accordance with another embodiment of the invention. The 
system of Fig. 7 includes a BTS 160 with a beamforming module 
166 and a plurality of antennas 168, 170, 172, 174, and two UEs 
162, 164, each having a pair of antennas 176/178, 182/184 and a 

15 decoder 180, 186. The beamforming module 166 includes 

respective beamformers 188, 190 and signal combiners 192/194, 
196/198 for each UE 162, 164. The beamformers 188, 190, and 
possibly the signal combiners 192, 194, 196, 198, are 
preferably implemented using signal processing software in 

20 conjunction with either a DSP or a general-purpose processor at 
the BTS 160. The decoders 180, 186 are also preferably 
software-based. 

Closed-Loop Transmit Diversity (CLTD) , with different 
configurations for the two UEs 162, 164, is provided by feeding 

25 back channel state information shown as H x , H 2 from the UEs 
162, 164 to the BTS 160. It should be appreciated that the 
combined 4x2 (M=4, N =2 , U =2) multi-user system of Fig. 7 is 
one illustrative example of a MIMO system to which the present 
invention may be applied, and that the invention is not limited 

30 thereto. This will become apparent from the following 
generalized analysis of multi-user MIMO with CLTD. 
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The BTS 160 has K ch spatial channels and M =4 
antennas 168, 170, 172, 174, with K ch < M . In a multi-user 
system, there exist only K chi channels to the i th UE . Signal 
vectors s(i) of symbols intended for the i th UE are multiplied in 
5 the beamformers 188, 190 by respective elements of 

F = [F il \F i2 \...F iU) ] , which is a spatial coding matrix of spatial 

coding weights. Each element F il) is a personal beamforming 
matrix for i th UE with MxK chi dimension and satisfying 

tr{F {i) F (i)T } = tr{F u)T F (i) } = P s , 2,.., U (29) 

10 where 

^{•} is the trace of a matrix. 

As shown, signals output from each beamformer 188, 
190 are combined in the signal combiners 192, 194, 196, 198 and 
output to respective antennas 168, 170, 172, 174. 

15 In this case, a constant and equal transmitted power 

is applied to signals for all UEs. The vector signal received 
by the i th UE is described by 

p(0 =H V) F 0)j(o + H U) ^ F {n) s M +77 (0 =H {i) Fs+fj (i) (30) where 

71=1, 

s u) =[s l 0) ,s 2 V) ...s Ktti 0} ] T is a K ch .-dimensioned vector of symbols 
20 of the i th UE; 

s = [s { ,s 2 ...ShY is a K ch -dimensioned vector of symbols of all 

UEs; 
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u 

K ch =^K chi is the total number of channels for the BTS; 
1=1 

J5 (0 = [y x {, \y 2 0) ...y N 0) ] T is an Af -dimensioned complex vector of 
signals received at the i th UE; 

H {l) is an NxM -dimensioned matrix of complex channel 
5 gains from the BTS to the i th user; 

7j U) =[?j l it \?j 2 0) ...rj N 0) ] T is an N -dimensioned complex vector of 
noise of observation for the i th user with zero mean and 

R i0 =2a TJ (l)2 I N covariance matrix. 

In multi-user systems in which UEs have different 
10 numbers of antennas, N would be replaced with N f above, where 
N g is the number of antennas at the i th UE. 

CLTD multi-user MIMO can be considered as the 
optimization of signal detection matrices. When the number of 
transmitting antennas is the same as the number of receiving 
antennas, so-called transmit and receive channel reciprocity 
exists. In this case, optimum weighting coefficients for 
beamforming for each UE are calculated at a BTS. The computed 
weighting coefficients are then used for radiation of a signal 
by the BTS. If the total number of receiving antennas of all 
UEs is equal to the number of transmitting antennas of the BTS, 
the full division of each transmitted signal, directly, on UE 
receiving antennas is satisfied. 

Now consider a reverse channel problem for a single 
user, to determine a CLTD matrix for data transmission from an 
25 i th UE, with N. antennas on K chi parallel channels. The virtual 



15 



20 
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reverse channel signal, observed at the BTS 166, is described 
by 

j)<o = ^(/)^(0|(0 + (31) 

where 

5 // (,) = [// (,) | is an A/ x Af. -dimensional channel matrix of the 

virtual reverse MIMO channel; 

F {0 is the optimal virtual CLTD space time coding matrix 
of N^K^ dimension; and 

the w/v " symbol indicates matrices and vectors associated 
10 with the virtual reverse channel. 

There are several ways to construct the F {,) matrix. 
Consider first a singular decomposition of the channel matrix 

H , H = UAV H , where U and V are unitary matrices with MxM 
and N i xN i dimensions, respectively, and A is a non-negative 
15 diagonal matrix with MxjV. dimension. The squares of diagonal 
elements of A are equal to eigenvalues of the HH' matrix. The 
columns of U are eigenvectors of the HH' matrix, and the 
columns of V are also eigenvectors of the HH' matrix. 

The optimum value of F can be shown to be 

20 F = V<S> , (32) 

where V is a matrix with N i *K chJ dimension, constructed from 
K ch . columns of V , and $ is a diagonal matrix with K ch i xK ch j 
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dimension having non-negative diagonal elements satisfying the 
following condition: 



(33) 



The diagonal elements of matrix O determine channel power 
allocation. A uniform power allocation gives 



h/=PJK ch4 . 



(34) 



Some possible alternative versions of power 
allocation include: 

MMSE criterion 



4>k,k = 2o V 



M 1 



(35) 



Minimum Symbol-Error-Rate (MSER) criterion 



2 2<V 



*k,k 



log 



( £ ^ 

hk y k 

K 2 < J 



; and 



(36) 



Maximum Capacity and Information Rate (MCIR) criterion, 
also commonly known as the water-filling algorithm 



2a 



2 V 



"k,k 



(37) 



where 
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(•) + =max(*,0) = l(|.| + #); 

4 kk = X kk 2 are eigenvalues of the HH f matrix, and X kk are 
diagonal elements of the A matrix; and 

ju is a factor that is selected to define each criteria. 

5 After the CLTD matrix is constructed at the 

transmitter, equation (31) becomes 

j> (0 = H F a) i (i) + fj {i) i = l 9 2 9 ...U , (38) 

where 

H F i0 = (H (i) F (i) )/^2a rj . 2 is a matrix of the virtual reverse 
10 MIMO channel with MxK chi dimension; and 

77=^(0,/^) is an M -dimensioned complex vector of a virtual 
noise observed at a BTS, with zero mean and R n =I M covariance 
matrix . 

A personal beamforming matrix F 0) can thereby be 
15 constructed in a closed loop fashion for each individual user 
in the absence of other users. However, in a multi-user 
scenario, the presence of inter-user MIMO interference prevents 
such a straightforward user-specific personal beamforming 
approach. An embodiment of the invention provides a solution 
20 for the multi-user CLTD by using the MMSE criterion to minimize 
the inter-user MIMO interference, and a network solution for 
optimizing the multi-user MIMO allocation. 
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To integrate signals of all users into a virtual 
model, it is possible to re-write the virtual reverse MIMO 
channel model for multiple users, as given below: 

? = H F i + fj , (39) 

where 



is a K ch -dimensioned vector of 



u 

symbols, and K ch .— ^ K chi ; and 



H F = \h f {1} 9 H f {2 \...H f ** J is an M x K ch -dimensioned matrix of the 
reverse virtual MIMO Channel. 

It should be noted that K ch < K ch < M , i.e. the number 
of estimated symbols is less than or equal to the number of 
received signals. In this situation, a very effective 
estimation can be carried out using, for example, a linear MMSE 
algorithm, as follows: 

§ = Gp , (40) 
where G = (H F H F +I, )' 1 H F = H F {H F H F ' + I k ) _1 



or 



5 (/) =G (,) J> , (41) 



where 
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G w = H F i0 (H F H F +I k )-' , i = l,2,. .U; and 



G = 



G 0) 



Using a principle of identity (or duality) of 
receiving and transmitting channels, an optimum demodulation 

5 matrix G is preferably used for generating F - [F (l) ,F (2) ,...F <£/) ] , 
the CLTD or beamforming matrix, at a BTS . In this case, 
personal beamforming matrices are preferably determined by 



F"=JP 



= JP. 



(42) 



m=l n=\ 



where g mn is (m,n) th element of G {) . 

10 The model of an observed (received) signal at the 

input of the UE of the i th user can be written as 



U 



= H F i0 s+fj (i) 



(43) 



where H F u ' n) = H (i) F (n) , H F W =[H F m \...H p vm '\ . In some communication 

systems, the BTS can measure the H F {, ' n) n = \,2...U matrices. In 
15 other systems, the UEs feed back channel matrices to the BTS. 
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The BTS can compute the personalized demodulation 
matrix for the i th user as 

L-i ( J f (44) 

and sends to each UE its respective demodulation matrix G (,) . 

5 Integrating all transformations for calculating the 

CLTD matrix (personal beamforming matrices F (i) at a BTS side 

and personal beamforming matrices G in at UE side) , yields the 
following algorithm to achieve CLTD based multi-user MIMO 
transmission in accordance with an embodiment of the invention: 

10 At a BTS side : 

g"=h/[h f h f ' +i Nl ]~ l 
h (i) = [h w ] 

At a UE side: 



6<o =H «Af j H««H^' +2* n 2 I Ni 

L«=i 

H V,n) = H (i) F (n) 



(46) 
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In the above multi-user system, the following 
constraints are preferably satisfied: 



1—1 



The first constraint of (47) specifies that the total 
5 number of parallel channels should not exceed the number of BTS 
antennas, which actually determines the number of parallel 
spatial channels. The second constraint specifies that the 
number of receiving antennas at a UE should be greater than or 
equal to the number of parallel spatial channels assigned to 
10 it. These conditions allow the use of linear methods to 

construct the personal beamforming matrices at both the BTS and 
UE for all users. 



For TDD (Time Division Duplexing) communications, all 
the computing can be done at a BTS. The BTS determines all 
15 relevant parameters, calculates both the BTS and UE beamforming 
matrices, and feeds back a personal beamforming receive matrix 

G (> to each UE. 

For the FDD (Frequency Division Duplexing) case, all 
the UEs can determine and feed back the initial SVD (Singular 

20 Value Decomposition) beamforming matrix F {i) to the BTS, which 

then jointly integrates all F li) matrices to compute the 

dedicated beamforming transmit matrix F U) for each UE . The BTS 
then computes and sends a respective beamforming receive matrix 

G ii} to each mobile terminal. 



25 Embodiments of the above CLTD-based multi-user MIMO 

exhibit performance advantages relative to conventional open- 
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loop solutions. For the purposes of comparison, the following 
simulation conditions were used: (1) R=l/2 turbo coding, block 
length 1280 bits, (2) QPSK (Quadrature Phase Shift Keying) 
modulation, (3) MMSE receiver for all schemes., and (4) ideal 
5 channel feedback . 

Before proceeding with a discussion of simulation 
results, each of the conventional technologies with which 
comparison is made will be briefly described. 

Null beamforming is an open-loop scheme where the 
10 number of transmitting antennas is equal to the total number of 

u 

receiving antennas, namely ^K chi =M r and K chi = N. , such as 

i=i 

shown in Fig. 8. Therefore, it is possible to perform null 
beamforming to reduce the inter-user interference of the BTS 
antenna emissions to the other users. The configuration of 
15 2x8x2 for a 4 user environment was considered in the 
simulations presented below. 



In one known open-loop multi-user BLAST technique, 
the number of transmitting antennas is the same as the number 
of receiving antennas for each UE, as shown in Fig. 9, and the 
20 number of channels is equal to the number of transmitting 
antennas. A UE receive signal can be expressed as 

y(0 _ / P * H i0 s +n (i) , where s is an M -dimensioned vector of 

v K ** 

symbols. At N t >M, it is possible to use an MMSE algorithm 
for demodulation of the entire s vector. Therefore, from the 
25 demodulated s vector, the symbols transmitted to a particular 
user are extracted. This scheme is generally known as Multi- 
user BLAST. We evaluate the configuration for such a system 
with dimension 2x8x8 and U = 4 . 
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Let K chi =M/U and Af. >M. In this case, it is 
possible for a BTS to group together the signals of each user 
into units with M symbols to thereby obtain units to 
sequentially transmit in a round robin fashion. At a UE side, 
5 the i th time slot is demodulated by the i th user with an MMSE 
decoder. This type of system in shown in Fig. 10. The 
performance of such a system is the same as a point-to-point 
BLAST system, with the compression of the duration of radiation 
of a signal of each user by a factor of l/U in time without 
10 increasing the transmit power. We evaluate the configuration 
2x8x8 and U = 4 for such a system. 

Figs. 11-13 show simulation results for embodiments 
of the invention and the above known technologies. 

Fig. 11 shows a plot of BLER versus SNR for an 
embodiment of the invention and several known communication 
schemes. It can be seen from Fig. 11 that for 8 transmit 
antennas at a BTS and 4 concurrent users each with 2 receive 
antennas at each UE, a CLTD scheme according to an embodiment 
of the invention has 19dB gain over open-loop null beamforming. 
This simulated CLTD scheme also achieves virtually the same 
performance as the 4 users open-loop multi-user BLAST each with 
8 receive antennas. In this case, a reduction of 6 receive 
antennas at each UE can be realized. Given 8 receive antennas 
at the UE side, the proposed CLTD scheme achieves 5dB gain over 
the conventional open-loop TDM BLAST, 9dB gain over multi-user 
BLAST, and 3 . 8dB gain over the closed-loop TDM BLAST. 

Fig. 12 is a plot of BLER versus Eb/No for an 
embodiment of the invention and several known communication 
schemes. The simulated CLTD scheme from which the plot of Fig. 
30 12 was generated is clearly superior to point-to-point MIMO 
(open-loop BLAST) and null beamforming. 



15 



20 



25 
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Fig. 13 is a plot of BLER versus SNR for several 
embodiments of the invention, and demonstrates the scalability 
of multi-user MIMO and MISO (Multiple Input Single Output) 
configurations according to embodiments of the invention. For 
5 the different configurations indicated in Fig. 13, all users 
achieve substantially the same level of QoS (Quality of 
Service) . 

Numerous modifications and variations of the present 
invention are possible in light of the above teachings. It is 
10 therefore to be understood that within the scope of the 

appended claims, the invention may be practised otherwise than 
as specifically described herein. 

Of course it is to be understood that in a given 
application, specific parameters may change. For example, 

15 different numbers of users, transmit antennas, and receive 
antennas may change the particular derivation details and 
equations above. However, adaptation of the above and further 
embodiments of the invention to other types and dimensions of 
systems than those explicitly described will be apparent to 

20 those skilled in the art from the foregoing. 

It should also be appreciated that references to 
transmitting or sending signals is not intended to limit the 
invention only to embodiments in which signals are transmitted 
exactly as generated, without any further processing. For 
25 example, signals may be compressed or otherwise processed prior 
to transmission, or stored for subsequent transmission at a 
later time. 



